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We present a new direct spectroscopic calibration for a fast estimation of the stellar metallicity [Fe/H], This calibrations was computed 
using a large sample of 451 solar-type stars for which we have precise spectroscopic parameters derived from high quality spectra. The 
new [Fe/H] calibration is based on weak Fe I lines, which are expected to be less dependent on surface gravity and microturbulence, 
and require only a pre-determination of the effective temperature. This temperature can be obtained using a previously presented 
line-ratio calibration. We also present a simple code that uses the calibrations and procedures presented in these works to obtain both 
the effective temperature and the [Fe/H] estimate. The code, written in C, is freely available for the community and may be used as 
an extension of the ARES code. We test these calibrations for 582 independent FGK stars. We show that the code can be used as a 
precise and fast indicator of the spectroscopic temperature and metallicity for dwarf FKG stars with effective temperatures ranging 
from 4500 K to 6500 K and with [Fe/H] ranging from -0.8 dex to 0.4 dex. 
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1. Introduction 

The derivation of stellar parameters is of extreme importance for 
several fields of astrophysics. Our knowledge of the fundamen- 
tal parameters of stars, such as the mass and radius, depends di- 
rectly on the precision that we can achieve when measuring the 
stellar atmosphere parameters. Parameters that we can directly 
derive from observations include the effective temperature, sur- 
face gravity, and metallicity. 

The two major techniques that are normally used to deter- 
mine these stellar parameters are photometry and spectroscopy. 
There are many calibrations to help us determine parame- 
ters based on the first technique. The IRFM method is one 
of these methods, that nowadays, is very commonly used be- 
cause of its reliability in determining th e effective temper- 
ature of a wide rang e of stellar t y pes ([B lackwell & Shallis, 
[!977t iBlackwell et all [l979l [l980l iBell & GustafssonTi l989; 
Gonzalez Hernandez & Bonifacio, 2009; Ramirez & Melendez, 
12005b lAlonso et all Il996t ICasagrande et all 120061 120101) . 
Other parameters such as the metal content can also be de- 
rived from photometric calibrations (iNordstrom et all [2004; 
ISchuster & Nisserl 1 19891) . These parameters can then be used 
to derive more fundamental parameters, such as mass and ra- 
dius u sing calibrations similar to those described in lTorres et all 
(2010). 

The determination of the spectroscopic parameters is not 
straightforward. A careful analysis is necessary for each stel- 
lar spectrum. This can be very time-consuming if you need pre- 
cise parameters and chemical abundances for several elements. 
One of the most difficult aspects of a spectroscopic analysis is, 



* Based on observations collected at La Silla Observatory, 
ESO, Chile, with the HARPS spectrograph at the 3.6-m tele- 
scope (072.C-0488(E)). The code is available for download at 
www.astro.up.pt/~sousasag/TMCalc 



undoubtably, the continuum determination. It can be very diffi- 
cult to derive correct position of the continuum, especially when 
dealing with poor quality spectra (considering data of low signal- 
to-noise (S/N) and/or low spectral resolution). Moreover, this 
difficulty in determining the continuum position also depends 
on the spectral type of a star. This is because cooler stars tend to 
have more spectral lines, increasing the amount of blending and 
leading to the "dilution" of the true continuum. The same hap- 
pens for different wavelength regions of a given stellar spectrum 
where the line density strongly increases as you move to shorter 
wavelengths. This problem can be strongly reduced using au- 
tomatic tools that are consistent in determining the continuum 
position, eliminating the errors/offsets, caused by a subjective 
determination, that are present when using interactive tools. 

Another problem that has an important impact on the deriva- 
tion of spectroscopic stellar parameters is the selection of the 
lines and their respective atomic parameters. There are two typ- 
ical choices made by spectroscopists. On the one hand, one can 
choose to use for each line the atomic parameters defined by a 
laboratory analysis. This method can have errors as large as 10- 
20 %. On the other hand, one can adopt new atomic -parameter 
values determined based on the observed lines in the spectra of 
a reference star (typically the Sun) and then assume its "well" 
known spectroscopic parameters. This second option allows us 
to perform a differential analyses using another star as a refer- 
ence, and, therefore, when the Sun is used as a standard candle, it 
is very suitable for solar-type stars (FGK). However, the method 
becomes imprecise for sufficiently cool and hot stars. 

This second option, which is usually referred to in the liter- 
ature as a differential spectroscopic analysis, was combined pre- 
viously with automatic codes su ch as ARES (ISousa et all 1 2007) 
applied in our previous studies (ISousa et all 120081120101 [201 1). 
This allowed us to derive in a systematic and homogeneous way 
spectroscopic stellar parameters for more than 1000 FGK stars. 
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These were t hen used to derive abundances for a large number 
of elements (iNeves et all l2009t lAdibekvan et all 1201 ll) . In all 
these works, the derived stellar parameters were proven to be 
compatible with others derived in different works using a range 
of independent methods. 

In this work, we make use onc e again of the data of 451 
stars presented in lSousa et all (120081) . which consist of very-high 
quality spectra of both high resolution and high S/N, to perform 
a new direct spectroscopic calibration of the stellar metallicity 
[Fe/H]. The calibration is based on weak Fe I lines that depend 
mainly on temperature and iron abundance and are less depen- 
dent on other parameters such as the surface gravity and micro- 
turbulence. Therefore, it only requires a pre-determination of the 
effective temperature. For the determination of the te mperature, 
we ca n use the line-ratio calibration presented in Sous a et al] 
(2010) and combine both calibrations to build a simple and fast 
code that allows a precise estimation of both the effective tem- 
perature and [Fe/H] . 

In Section 2 we recall and discuss the main features of the 
line-ratio calibration code, presenting some tests done in previ- 
ous studies to demonstrate its consistency in the particular pa- 
rameter space of the calibration. In Section 3, we define the new 
[Fe/H] spectroscopic calibration and explain how we derive it. 
We also present a simple procedure to derive the final value 
of [Fe/H] obtained from all the individual line calibrations. In 
Section 4, we show some simple tests performed on both the 
calibration sample and an independent large sample of solar-type 
stars. In the final Section 5, we summarize this work. 

2. Temperature based on a line-ratio calibration 

The code presented in this work is inspired by a previously pub- 
lished effective-temperature calibration bas ed on line ratios o f 
several spectral lines of different elements (ISousa et al.. 2010). 
This calibration was presented as an excellent tool for deter- 
mining automatically and quickly a spectroscopic effective tem- 
perature, and can be easily used to confirm a spectroscopic 
temperatu re determined using the "standard" procedures (see 
alsolGrav & Johansonl[T99llGrayL[T994l;lGrav & Brownll200lt 
lGravLl2004UKovtvukh et alll2003l) . This was presented as a pos- 
sible extension of the ARES code, which can automatically mea- 
sure t he equivalent wid t hs (EWs) of weak ab sorption spectral 
lines (ISousa et alll2007l) . lDesidera et ai1(l201 ll) showed that this 
line-ratio calibration can return good estimates of the tempera- 
ture, even for solar-type stars up to relatively high rotation rates 
(up to 18 Km/s). This excedes the abilities of typical EW meth- 
ods owing to the increase in the number of blends with higher 
rotation of the stars. 

In Figure [1] we com pare the temperatu res derived with the 
line-ratio calibration of ISousa et al l (1201 Oh with those derived 
using the standard spectroscopic procedure (ISousa et al. L I2008L 
1201 ll) . The top plot shows this comparison for the sample of stars 
used to calibrate the line ratios. As expected, the consistency 
is very good. This plot is important to illustrating the typical 
dispersion that we obtain using the line-ratio calibration. 

In the bottom plot of the same figure, we show the same 
comparison for an "independent" sample of stars (independent 
in the sense that these stars were not used to compute the cali- 
bration of each line ratio). It is clear that the comparison is also 
consistent. Both plots were presented in previous works, and 
they are presented here to recall that the temperature inferred 
from the line ratio is consistent with our standard spectroscopic 
method, within the ranges defined for the calibration. Therefore, 
we can use the line-ratio calibration in order to determine the 
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Fig. 1. In the top panel, we show the comparison between the 
effective temperatures derived from our standard spectroscopy 
analysis and the ones derived using the line-ratio code for the 
sample used to deri ve the calibration ( labelled Cal. Sample in 
the plot adapted from lSousa et al.l(l2010l) ). In the bottom plot, we 
sh ow a similar plot b ut this time for the sample of stars presented 
in lSousa et al.l d201 ll) . This corresponds to a sample that is inde- 
pendent of the derived calibration. In both panels, we present the 
mean difference and the standard deviation for the comparisons. 



temperature of a star, and then examine the strength of the iron 
absorption-lines and quickly find a calibration of [Fe/H]. 

3. [Fe/H] calibration 

3. 1 . Selection of lines 

To obtain a calibration that help us determine [Fe/H], we need to 
compile a list of iron lines that are observed for solar-type stars. 
For this step, we s tarted with the "stable" line-list presented in 
~aTT (l2008l) . which is a list of iron lines chosen to ensure 



ISousa et" 

that they are appropriate for an automatic determination of EWs 
with ARES. The next step was to select lines that should be in- 
dependent of both surface gravity and microturbulence. In other 
words we needed at this point, to select lines that depend mostly 
on temperature and [Fe/H]. 

We can easily ensure the independence of the calibration of 
the surface gravity by selecting only the iron lines that corre- 
sponds to the neutral state (Fe I). This is because the ionized 
iron lines (e.g. Fe II) depend strongly on the surface gravity, in 
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Fig. 2. The sequence of plot from top to bottom shows the re- 
moval of outliers for the iron line at 6820.37 A. 



contrast to the neutral iron lines (lGravlll992l) . We therefore take 
out all the Fe II lines from the original list. 

Finally, the microturbulence parameter may have a strong 
impact on the strength of a spectral line, particulary for stronger 
lines. Weak lines are known to be reasonably independent of 
microturbulence. 

One problem that we face in defining our linelist is that the 
strength of each individual line depends strongly on the effective 
temperature and [Fe/H]. We can therefore only eliminate a line 
after its line strengh has been measured and checked to be in the 
regime where it is independent of microturbulence. Therefore, 
for this work we consider our cut at 75 mA. This value seems to 
be reasonable taking into account not only the microturbulence 
dependence, but also considering that lower values will strongly 
reduce the number of lines acceptable for the calibrations. 

We also wish to note that there is a kno wn dependence of 
the microturbulence on the temperature (e.g. Pilachowski et al.L 
1996: iFor & SnedenLl2010l) . Therefore, the microturbulence de- 
pendence seen for the stronger lines may be removed signifi- 
cantly by the temperature fitting in our calibrations. Together 
with the restrictions that we impose on the line selection we 
can strongly eliminate any microturbulance dependence from 
our calibrations. 

3.2. The [Fe/H] calibration 

To derive the calibration for [Fe/H], we again used the sample of 
stars presented before for the temperature line-ratio calibration. 



This sample is composed of 451 solar-type stars that were all 
observed with HARPS to obtain high-resolution spectra (R ~ 
1 10000). Thes e stars were part o f the HARPS survey for extra- 
solar planets (Mavo ret all 1201 it) . The S/N for this sample varies 
from 70 to 2000, with 90% of the spectra having a S/N higher 
than 200. It is a very well-established sample with homogeneous 
spectroscopic parameters that have been compared with other 
independent methods. 

To ensure consistency with the previous line-ratio calibra- 
tion, we used the same EW measurements. 

For each line, a calibration was computed following the 
equation: 



EW =C + Ci * [Fe/H] + C 2 * T eff + 
C 3 * [Fe/H] 2 + C 4 * T 2 eff + 
C 5 * [Fe/H] * T eff . 



(1) 



We choose to use this equation in order to have an easy (and 
quickly attainable) representation of the physical dependence 
between the line strength, the effective temperature, and the iron 
abundance. To extract the iron abundance from this equation, we 
perform a simple invertion of [Fe/H]. 

In Fig. [2] we show the results of an iterative process for the 
determination of the calibration coefficients for the iron line at 
6820.37 A. In this process, we first select the lines (stars) ac- 
cordingly to the discussion in Sect. 3.1, i.e. removing all points 
whith EWs larger than 75 mA (each point is the EW for a star 
in the sample). At this step, we also choose to remove the lines 
with EWs smaller than 20 mA. This is because for these small 
lines we have a larger relative error in the EW. We then fit the re- 
maining points with the equation and perform the first outlier re- 
moval, choosing a 2 cr threshold. We repeat the process one more 
time to eliminate extra outliers. These outliers are mainly due to 
poor automatic EW measurements that are typically related to 
either a bad continuum determination or strongly blended lines. 
The final fitted coefficients are kept, together with the final num- 
ber of stars used, the slope of the direct comparison between the 
calibrated and the spectroscopic [Fe/H], and the respective stan- 
dard deviation of the fit. 

This procedure was performed for all the lines in the initial 
iron line-list. Finally, we made a careful selection of the lines, 
considering the results of each individual fit. This selection was 
performed to ensure that each calibration satisfies the following 
conditions: 

1 . The direct comparison between the calibrated and the spec- 
troscopic [Fe/H] has a slope within 3% of the identity line. 

2. The standard deviation of the individual line calibration is 
less than 0.04 dex. 

The values used in these two conditions were chosen to keep 
a significantly large number of lines in order to increase the sta- 
tistical meaning of the final metallicity derivation. In this case 
we want to note we choose not to neglect lines that have a given 
minimum of the final number of stars used for the individual 
fit. The reason for this is that we remove lines of the fit ac- 
cordingly to their strength, which means that depending on the 
temperature and metallicity of a star, a given line will be either 
stronger or weaker in specific space regions of these two param- 
eters. Therefore, the only concern that we have in these cases is 
to ensure that we take into account that each individual calibra- 
tion is only valid within each individual parameter space. 
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Table [2] lists a sample of the total of 149 lines that have 
passed this process and can be used for our [Fe/H] determina- 
tion. The full table is available in its electronic format. 

3.3. [Fe/H] estimation from the line calibrations 

From the calibrated lines, we can now derive a final value for the 
global [Fe/H] for a given star. The procedure that we propose 
here is the same as that p resented for the tem perature determina- 
tion based on line ratios ( ISousa et"ail 12010). We summarize this 
procedure in the following items: 

- First we compute the [Fe/H] determination using each cali- 
brated line from Table [2] 

- Secondly, we select the calibrated lines accordingly to the 
limits of each individual calibration. In this step, we choose 
to increase the limits by 100K in both directions. The errors 
coming from the temperature line-ratio calibration are of this 
order of magnitude and therefore we wish to guarantee that 
we do not remove lines that can give a reasonable estimation 
of [Fe/H]. 

- Finally, we compute the weighted average of the [Fe/H] re- 
sults, considering the standard deviation of each individual 
calibration. 

This procedure is repeated twice with a 2 cr outlier removal, 
eliminating in this process the EWs that were not (for any rea- 
son) performed correctly. 

3.4. Errors 

The easiest error estimate that we can extract from our procedure 
is to assume the dispersion in the values given by all the individ- 
ual calibrations. If one considers that each individual calibration 
is independent of all the others, we can divide the dispersion 
by the square root of the final number of individual calibrations 
used. 

It is wise, however, to include the error in the temperature. 
To do this, we choose to use a straightforward estimation of this 
error, which is to derive the error in [Fe/H] using the limits given 
by the temperature error (1-cr). The final error is obtained by 
evaluating the quadratic sum of the two sources of errors. 

3.5. Using TMCalc to estimate temperature and [Fe/H] for 
solar-type stars 

The calibration presented here is only useful when you have a 
temperature estimation. We therefore developed a free code, im- 
plemented in C language, that combines both the line ratio cal- 
ibration and the [Fe/H] calibration presented here. The code is 
available online at the ARES web pag^J The code comes with 
a simple shell script "TMCalc", which can be used as the driver 
to run the C code with an easy shell-command line. This makes 
the code ready to be included in any kind of a spectral analysis 
pipeline. The only requirement is to have the spectrum (e.g. com- 
ing from the pipeline) in a format ready to be used with ARES. 

4. Testing the code 

4.1. The calibration sample 

Figure [3] shows the results obtained when applying the code to 
the calibration sample. Here we compare the calibrated [Fe/H] 

1 http://www.astro.up.pt/~sousasag/ares) 




Fig. 3. Direct comparison between the calibrated [Fe/H] and the 
spectroscopic [Fe/H] for the sample used to compute the calibra- 
tions. 



against the spectroscopic [Fe/H]. The result is consistent, thus 
can be expected since in this case the sample is the same as that 
used to compute the line calibrations. 

In Fig. [4] we try to identify any dependences of the value of 
[Fe/H] derived from the calibration on other spectroscopic pa- 
rameters such as temperature, surface gravity, and spectroscopic 
[Fe/H]. It is possible to see from this figure the ranges of spec- 
troscopic parameters for which this calibration is valid. It is clear 
that these come directly from the star properties in the sample, 
which is typically composed of solar-type stars with effective 
temperatures ranging between 4556 K and 6403 K, and with sur- 
face gravities typical of dwarf stars and a few sub-giants ([3.68, 
4.62]). All these stars have metallicities of around solar, ranging 
from -0.84 dex to 0.39 dex. 

It is clear from Fig.[4]that there is no clear trend, except for 
the spectroscopic [Fe/H] itself. An offset is clearly visible for the 
metal-rich stars, but the maximum differences remain within 0. 1 
dex, which can nevertheless be significant. This offset is a result 
of some extra dependence that the individual line calibrations 
were unable to extract, or is merely an artifact of the fit, since 
the offset is only observed in the upper limit to [Fe/H]. Instead 
of applying a higher order polynomial for each individual cali- 
bration and trying to eliminate this offset, we choose to perform 
a simpler and more direct correction. From the figure, we can see 
that we can fit the points with a simple horizontal line for [Fe/H] 
< -0.2 (dashed line). This corresponds to a calibrated [Fe/H] 

0.187, given that the horizontal line has a value of ~ 0.013. 

For higher [Fe/H] values, we use a second-order polynomial to 
consistently fit the offset (dashed-point line). We therefore use 
the following equations to perform the correction for the cali- 
brated [Fe/H] in these two regimes: 

[Fe/H] cor =[Fe/H]cai- (0.013) 

if ([Fe/H] cal < (-0.187)), 

[Fe/H] cor =[Fe/H] ca , - (-0.000726960 (2) 

- 0.0950966 * [Fe/H] ca , 

- 0.128329 *[Fe/H] 2 cal 

if ([Fe/H]cai> (-0.187)). 
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Fig. 4. Difference between the calibrated [Fe/H] and the spec- 
troscopic [Fe/H] and its dependence on the other spectroscopic 
stellar parameters such as the effective temperature (top panel), 
the surface gravity (middle panel), and the spectroscopic [Fe/H] 
(bottom panel). The dashed and dashed pointed lines are fits to 
the date described in detail in the text. 



This equation is applied at the end of the procedure to correct 
the offset that can be seen in Figure [4] 

In Figure[5] we show the effect of the correction for the sam- 
ple used for the calibration. The dispersion in slightly smaller 
and the offset is close to zero as expected. However, this is still 
the sample of stars used to derive the calibration, hence the con- 
sistency of this result is expected. What is interesting to show 
at this stage is the "small" dispersion in the evaluated [Fe/H], 
which is typically around ~ 0.03 dex. This value can be used as 
a reference to indicate the quality of the final corrected calibra- 
tion. 

4.2. Constrains on line strenghts 

As discussed before, we choose lines of a specific range of 
strengths: weak lines in order to avoid any dependence on the 
microturbulence, but not those that are weak to avoid the intrin- 
sic errors in the measurements of very weak lines. The reader 
may be concerned about these restrictions and the possible sys- 
tematics that they can produce on the described calibrations. For 
instance, for low temperature stars, these restrictions will tend to 
eliminate the stronger lines in the metal-rich stars, while in hotter 
stars, different lines (that are weaker) will be eliminated for the 
metal-poor stars owing to the poorer constraints. Even without 
our poorer constraints, the iron lines tend to disappear from the 
spectrum in these cases, since they are metal-poor. Therefore, 
for different types of stars, we use different sets of lines, which 
could introduce systematic errors into the calibrations. 
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Fig. 5. Direct comparison between the corrected calibrated 
[Fe/H] and the spectroscopic [Fe/H] for the sample used to com- 
pute the calibrations. 

Figure [6] tries to identify any of the described possible sys- 
tematic errors. In the top panels of this figure, we show the 
same kind of direct comparison between the corrected calibrated 
[Fe/H] and the spectroscopic [Fe/H], but divide the sample into 
different temperature regimes. On the left we see the lower tem- 
perature stars (T e ff < 5000 K), and on the right the hotter stars 
(T e /f > 5800 K). The comparisons remains consistent with a 
small dispersion and no visible trend for the different temper- 
ature regimes. In the bottom panels, we present the number of 
individual calibrations used in each star. In these plots the de- 
crease in the number of lines used is clear for the metal-poor 
stars in both temperature regimes. For the cooler stars, the re- 
duction in the number of lines is also clear for the more metallic 
stars. This is the effect that we described before and is a result 
of neglecting the stronger lines to eliminate the microturbulence 
dependence. Although the number of lines is smaller owing to 
the discussed systematics, it is clear from the comparison that 
the results are consistent, and the selection of the lines using the 
indicated constraints on the line strength seem to produce good 
results. 

4.3. Large independent sample for direct comparison 

We only showed the test on a very well-defined sample, com- 
posed of stars with data of high resolution and high S/N, which 
was used to define the calibration. At this point, we wish to show 
an independent test, in the sense that we use a different sample 
of stars that was not used to compute the calibration. This "in- 
dependent" sample is composed of a total of 582 FGK stars for 
which the parameters were determined using the same homoge- 
neous spectros copic method. The details of these determinations 
can be seen in Sou sa et al ] (I2011h . The main difference regards 
the spectral quality of this sample: this sample typically has a 
lower S/N than the calibration sample. 75% of the stars have S/N 
lower than 200, while 90% of the stars in the calibration sample 
have data of S/N higher than 200. 

In Figure [7] we compare the standard spectroscopic deriva- 
tion result with the calibration presented here for this "indepen- 
dent" large sample of stars. 

The result of this direct comparison is very good. The mean 
difference is close to zero, indicating that there is no clear offset 
in the determined [Fe/H]. However, the dispersion (~0.04 dex) 
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Fig. 6. Direct comparison between the corrected calibrated [Fe/H] and the spectroscopic [Fe/H] for the sample used to compute the 
calibrations for two different temperature regimes: cool stars (Top left panel) and hotter stars (top right panel). In the bottom panel, 
we present the number of individual calibrations used in each star. 




Fig. 7. Direct comparison of the corrected calibrated [Fe/H] and 
the spectroscopic [Fe/H] for an independent sample of solar-type 
stars. The three grey squares indicate stars for which the code 
obtained results outside the calibrated [Fe/H]. 

is larger than for the previous comparisons. This is expected be- 
cause this sample has typically a lower S/N. In addition, this 
large "independent" sample has a few stars that lie outside the 
limits of the calibration. Amoung the 582 stars that belong to 
this sample, the code was able to derive values for a total of 
556 stars. Amoung these stars, there are still a few that are out- 
side the limits of the calibration. This can be clearly seen in the 
figure at lower metallicities where we have at least three stars 
with metallicities lower than -0.8 dex (plotted as grey squares in 
the Figure). The error bars for these three stars are large com- 
pared to those for the rest of the stars. This is due to the far 
smaller number of iron lines (about ten) used in the estimation 
of the calibrated [Fe/H]. The number of lines and their respec- 



tive errors derived from the procedure previously described can 
be used to reliably check whether the star is within the limits of 
the parameter space of the calibration. 

After this has been noted, the remaining results are very con- 
sistent and reveal that this calibration, both in temperature and 
[Fe/H], is precise enough to be used within its limits, to effi- 
ciently derive these two atmosphere parameters, either as a first 
approximation or as a check of other standard spectroscopic de- 
terminations. 

4.4. Spectral resolution and S/N 

To evaluate how this procedure works for different spectral res- 
olutions and different values of S/N, we computed several solar 
spectra, first with different levels of S/N and then with different 
levels of resolution. The original spectrum was obtained using 
the Kurucz Solar Atlas. The noise was created using a Gaussian 
distribution and we introduced the artificial inst rumental reso- 
lution using the "rotin3" routine in SYNSPEC0(Hubeny et al., 
1994). We first created ten different S/N levels to be included 
in the original spectrum (10, 25, 50, 75, 100, 150, 200, 300, 
400, and 500). We then created ten addicional different spectra 
with resolutions ranging from 10000 to 100000, corresponding 
to FWHM of 0.550 A and 0.055 A at the central wavelength 
5500 A. For this spectra, we included the noise to simulate spec- 
tra with S/N ~ 200, which more closely represent the real spec- 
tra. We generated the spectra in the region [4000 A - 7000 A ] in 
order to cover our linelist. 

For the spectra of different resolutions, the EW measure- 
ments were all performed with ARES using the same input pa- 
rameters (smoothder=4, space=3, rejt=0.994, lineresol=0.07), 
except for the spectrum with the lowest resolution for which we 
adapted a "smoothder" parameter of 12. For the different S/N 



2 http://tlusty.gsfc.nasa.gov/index.html 
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Fig. 8. Temperature and [Fe/H] determined using our calibra- 
tions for the simulated solar spectrum with different sets of res- 
olutions and S/N. We show two different error estimations, the 
thin black error bars representing only the dispersion in the in- 
dividual calibrations, and the thick grey error bar where we also 
take into account the number of individual calibrations used. See 
text for details. 



spectra, the ARES parameter "rejt" changed accordingly to the 
S/N of the spectra , using the recommended values discussed in 
Sousa et al] d201 lb . The results can be seen in Figure[8] The two 
top panels shows the results for the solar spectra of different S/N, 
the two bottom plots for the solar spectra of different resolution. 
We also present two of the errors estimations given by the code. 
The smaller and thicker grey error-bars represent the errors ob- 
tained by considering the dispersion in the individual calibration 
and the number of independent lines used. The thin black error 
bars presented in this figure are the ones obtained using only the 
dispersion in the individual calibrations, ignoring the number of 
calibrations used. The results are consistent even for the lower 
resolutions, where the strong increase in the errors originating 
directly from dispersion in the independent calibrations can be 
clearly seen. This proves that the calibrations presented in this 
work can be safely used for lower resolution spectra. In terms 
of the S/N, the results are also quite consistent down to S/N ~ 
25. For this low S/N, the temperature estimation starts to deviate 
from the expected one, although the errors are large and greater 



than the offset. We note that for the [Fe/H] even for very low S/N 
spectra the values are very close to the expected [Fe/H]=0. 

5. Summary 

We have presented a new spectroscopic line calibration to effi- 
ciently estimate the stellar [Fe/H]. This calibration depends on 
both the temperature and the strength of the iron lines. We make 
available a free code that combined with ARES and a previous 
temperature line-ratio calibration allows us to estimate both the 
spectroscopic stellar effective-temperature and the [Fe/H]. We 
tested this code with a large sample of solar-type stars and con- 
firmed that these calibrations are consistent within the parameter 
space defined for the calibrations. This code can easily be ap- 
plied to a spectroscopic data-analysis pipeline in order to quickly 
obtain precise estimations of these important spectroscopic pa- 
rameters. 

These calibrations should not replace the more precise 
standard spectroscopic analyses that are typically more time- 
consuming. These standard methods should definitely be used 
if one whishes to study stars individually. However, the tool that 
we present in this work can be very useful for determining pre- 
cise parameters for large amounts of data. Typical examples are 
the data sets of spectroscopic surveys, which are very large and 
for which people normally search for statistical trends among the 
propreties of the many survey targets. 
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5679.03 


1 


396e+02 


1.971e+01 


-3.645e-03 


-4.596e+0() 


-1.765e-06 


5 


075e-03 


4749 


6374 


-0.84 


0.33 
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Table 2. continued. 



Sousa, S. G. et al.: TMCalc - a fast code to derive Teff and [Fe/H] for FGK stars 



id 


stddev 


ns 


A 


co 




C2 


C3 


c 4 


C5 


iff 


eff 


fehmin 


fehmax t 


75 


0.017 


181 


5715.09 


1.853e+02 


2.684e+01 


-7.341e-03 


4.650e+0() 


-2.123e-06 


5.455e-03 


4649 


6403 


-0.83 


0.21 


76 


0.016 


350 


5731.77 


5.208e+01 


7.093e+00 


2.631e-02 


-1.336e+00 


-4.375e-06 


7.272e-03 


4649 


6403 


-0.84 


0.33 


77 


0.025 


123 


5738.24 


-1.048e+02 


1.380e+02 


6.458e-02 


2.746e+01 


-7.641e-06 


-1.874e-02 


4556 


5914 


-0.20 


0.39 


78 


0.020 


362 


5741.85 


-2.092e+01 


4.235e+01 


4.233e-02 


6.322e+00 


-5.725e-06 


3.309e-04 


4578 


6361 


-0.70 


0.39 


79 


0.022 


378 


5752.04 


-8.798e+01 


-1.295e+01 


7.141e-02 


-2.124e+00 


-8.080e-06 


1.065e-02 


4730 


6374 


-0.84 


0.35 


80 


0.018 


327 


5775.08 


-2.884e+()l 


-2.385e+00 


5.526e-02 


-2.22()e+0() 


-6.923e-06 


8.942e-03 


4728 


6374 


-0.84 


0.33 


81 


0.018 


287 


5778.46 


1.890e+01 


4.506e+01 


3.776e-02 


7.088e+00 


-6.442e-06 


-2.017e-04 


4556 


6161 


-0.70 


0.39 


82 


0.030 


387 


5793.92 


5.060e+01 


5.856e+01 


1.629e-02 


-2.570e-01 


-3.303e-06 


-3.114e-03 


4649 


6374 


-0.72 


0.39 


83 


0.031 


371 


5806.73 


-2.200e+02 


2.733e+01 


1.206e-01 


1.304e+00 


-1.263e-05 


4.599e-03 


4578 


6374 


-0.84 


0.34 


84 


0.032 


356 


5809.22 


2.561e+02 


1.174e+02 


-3.839e-02 


7.787e+00 


4.498e-07 


-1.101e-02 


4578 


6374 


-0.84 


0.35 


85 


0.014 


271 


5814.81 


-5.541e+01 


5.903e+01 


4.968e-02 


1.129e+01 


-6.253e-06 


-3.345e-03 


4556 


6227 


-0.54 


0.39 


86 


0.029 


99 


5827.88 


2.484e+02 


1.757e+02 


-6.769e-02 


2.119e+()l 


4.691e-06 


-2.744e-02 


4595 


5886 


-0.32 


0.39 


87 


0.021 


381 


5852.22 


5.740e+01 


6.228e+01 


1.907e-02 


5.890e+00 


-3.794e-06 


-2.186e-03 


4624 


6403 


-0.78 


0.39 


88 


0.015 


266 


5855.08 


-1.581e+02 


4.554e+01 


8.292e-02 


1.032e+01 


-8.931e-06 


-1.120e-03 


4578 


6227 


-0.38 


0.39 


89 


0.017 


386 


5856.09 


-1.141e+02 


2.988e+01 


7.359e-02 


6.299e+00 


-8.302e-06 


2.511e-03 


4624 


6403 


-0.70 


0.39 


90 


0.035 


95 


5902.48 


-2.58()e+02 


4.483e+01 


1.138e-01 


2.863e+01 


-1.150e-05 


-4.066e-03 


4723 


5914 


-0.07 


0.39 


91 


0.032 


319 


5916.26 


1.395e+02 


-8.026e+0() 


5.336e-03 


-4.553e+0() 


-3.450e-06 


9.420e-03 


4649 


6374 


-0.84 


0.33 


92 


0.030 


383 


5927.79 


-8.166e+()l 


-7.428e+00 


6.484e-02 


-2.926e+00 


-7.472e-06 


9.326e-03 


4578 


6403 


-0.78 


0.39 


93 


0.033 


393 


5929.68 


-6.996e+01 


-5.744e-()2 


5.875e-02 


-6.093e+00 


-6.873e-06 


7.4()2e-03 


4649 


6403 


-0.78 


0.39 


94 


0.022 


136 


5934.66 


9.081e+01 


-1.922e+01 


2.423e-02 


-4.132e+00 


-4.644e-06 


1.236e-02 


5023 


6374 


-0.84 


0.21 


95 


0.023 


292 


5956.70 


1.979e+02 


7.238e+00 


-4.371e-03 


3.304e+00 


-3.624e-06 


6.958e-03 


4820 


6403 


-0.84 


0.35 


96 


0.019 


305 


6027.06 


-3.348e+01 


-1.738e+01 


5.513e-02 


-1.116e+00 


-6.6()0e-06 


1.110e-02 


4671 


6374 


-0.84 


0.30 


97 


0.023 


190 


6056.01 


-1.434e+02 


-5.097e+01 


9.439e-02 


-3.554e+00 


-9.882e-06 


1.737e-02 


4578 


6374 


-0.84 


0.21 


98 


0.016 


142 


6078.49 


2.758e+02 


6.500e+01 


-4.065e-02 


3.100e+00 


1.085e-06 


-1.103e-03 


4879 


6374 


-0.84 


0.18 


99 


0.020 


390 


6079.01 


-7.753e+01 


3.215e+01 


6.533e-02 


9.778e-()l 


-7.603e-06 


2.745e-03 


4671 


6403 


-0.78 


0.39 


100 


0.017 


374 


6082.72 


1.050e+02 


4.141e+01 


1.286e-02 


5.201e+00 


-4.321e-06 


7.475e-04 


4624 


6361 


-0.84 


0.39 


101 


0.033 


375 


6089.57 


-4.388e+()l 


2.996e+01 


5.058e-02 


7.216e+00 


-6.371e-06 


2.512e-03 


4556 


6361 


-0.72 


0.39 


102 


0.020 


220 


6094.38 


-3.242e+()2 


3.546e-()l 


1.431e-01 


7.829e+00 


-1.445e-05 


7.125e-03 


4556 


6161 


-0.28 


0.39 


103 


0.017 


375 


6096.67 


3.028e+01 


6.386e+01 


3.210e-02 


3.028e+00 


-5.310e-06 


-2.434e-03 


4671 


6361 


-0.70 


0.39 


104 


0.026 


165 


6098.25 


-1.722e+02 


4.669e+01 


8.482e-02 


2.465e+01 


-9.018e-06 


-3.001e-03 


4556 


6161 


-0.25 


0.39 


105 


0.017 


358 


6151.62 


8.572e+01 


2.699e+00 


2.460e-02 


9.152e-01 


-5.321e-06 


7.643e-03 


4578 


6374 


-0.84 


0.35 


106 


0.031 


325 


6157.73 


-6.641e+01 


-2.480e+01 


6.653e-02 


-8.756e+00 


-7.682e-06 


1.155e-02 


4578 


6374 


-0.84 


0.35 


107 


0.036 


389 


6165.36 


-1.978e+02 


-3.319e+01 


1.046e-01 


-3.580e+00 


-1.081e-05 


1.340e-02 


4556 


6374 


-0.84 


0.39 


108 


0.015 


192 


6173.34 


2.765e+02 


2.456e+00 


-3.815e-02 


-1.829e+00 


3.782e-07 


7.731e-03 


4879 


6361 


-0.83 


0.30 


109 


0.016 


369 


6187.99 


1.513e+02 


6.428e+01 


-4.207e-03 


3.402e+00 


-2.372e-06 


-1.913e-03 


4750 


6403 


-0.84 


0.35 


110 


0.015 


159 


6200.32 


2.515e+02 


-1.202e+01 


-2.843e-02 


-2.718e+00 


-4.263e-07 


1.051e-02 


4900 


6361 


-0.83 


0.21 


111 


0.020 


239 


6220.79 


9.946e+01 


8.621e+01 


-2.926e-03 


1.944e+()l 


-1.916e-06 


-8.336e-03 


4556 


6161 


-0.56 


0.39 


112 


0.018 


336 


6226.74 


-4.656e+01 


6.611e+01 


5.3()8e-02 


1.105e+01 


-6.895e-06 


-3.410e-03 


4624 


6276 


-0.56 


0.39 


113 


0.019 


375 


6229.24 


-1.655e+00 


6.564e+00 


4.449e-02 


5.706e+00 


-6.467e-06 


7.753e-03 


4556 


6374 


-0.72 


0.39 


114 


0.038 


322 


6238.39 


-3.487e+02 


-3.509e+01 


1.158e-01 


1.462e+01 


-8.241e-06 


1.371e-02 


4831 


6361 


-0.67 


0.39 


115 


0.016 


360 


6240.65 


6.372e+01 


3.407e+00 


3.232e-02 


1.895e+00 


-6.022e-06 


7.794e-03 


4624 


6403 


-0.78 


0.35 


116 


0.016 


375 


6270.23 


2.415e+01 


1.398e+01 


4.272e-02 


2.571e+00 


-6.554e-06 


5.993e-03 


4624 


6374 


-0.84 


0.35 


117 


0.028 


391 


6315.81 


-1.179e+02 


1.939e+01 


7.977e-02 


3.463e+00 


-9.014e-06 


4.853e-03 


4578 


6403 


-0.70 


0.35 


118 


0.017 


140 


6322.69 


1.652e+02 


-8.515e+00 


1.986e-03 


-3.437e+0() 


-3.022e-06 


9.911e-03 


5023 


6374 


-0.84 


0.19 


119 


0.010 


94 


6358.68 


1.589e+02 


4.306e+01 


3.146e-02 


2.037e+00 


-7.674e-06 


4.734e-03 


5417 


6374 


-0.78 


0.21 


120 


0.036 


397 


6380.75 


-2.942e+()2 


-2.097e+01 


1.426e-01 


-5.316e+0() 


-1.429e-05 


1.144e-02 


4556 


6374 


-0.84 


0.35 


121 


0.018 


240 


6392.54 


-1.515e+02 


1.361e+01 


9.765e-02 


3.821e+00 


-1.188e-05 


4.910e-03 


4578 


5989 


-0.72 


0.39 


122 


0.024 


233 


6481.88 


-2.832e+01 


-4.473e+01 


6.769e-02 


2.801e-()l 


-8.946e-06 


1.604e-02 


4649 


6403 


-0.83 


0.30 


123 


0.027 


331 


6498.94 


2.621e+02 


-2.801e+01 


-3.054e-02 


-7.395e+0() 


-1.192e-06 


1.237e-02 


4649 


6374 


-0.84 


0.39 


124 


0.021 


247 


6608.03 


-1.104e+02 


3.988e+01 


8.333e-02 


5.617e+00 


-1.062e-05 


1.996e-04 


4578 


5989 


-0.72 


0.39 


125 


0.015 


219 


6609.12 


1.566e+02 


1.759e+01 


4.808e-03 


4.456e+00 


-3.536e-06 


6.485e-03 


4649 


6374 


-0.84 


0.30 


126 


0.023 


229 


6625.02 


1.285e+02 


6.232e+01 


6.883e-03 


8.233e+00 


-4.658e-06 


-3.701e-03 


4556 


5914 


-0.72 


0.39 


127 


0.018 


326 


6627.55 


-1.703e+02 


8.829e+01 


9.358e-02 


1.550e+01 


-1.023e-05 


-7.453e-03 


4578 


6260 


-0.50 


0.39 


128 


0.016 


380 


6703.57 


2.267e+01 


2.445e+01 


3.990e-02 


3.450e+00 


-6.453e-06 


4.010e-03 


4556 


6361 


-0.78 


0.39 


129 


0.018 


390 


6705.11 


-8.592e+01 


5.960e+01 


7.263e-02 


6.213e+00 


-8.585e-06 


-1.000e-03 


4556 


6361 


-0.78 


0.35 


130 


0.022 


229 


6710.32 


7.043e+01 


7.075e+01 


2.311e-02 


1.057e+01 


-5.682e-06 


-4.916e-03 


4556 


5917 


-0.72 


0.39 


131 


0.021 


292 


6713.05 


-2.097e+()2 


6.184e+()l 


1.048e-01 


9.106e+00 


-l.llle-05 


-3.766e-03 


4671 


6161 


-0.41 


0.39 


132 


0.017 


236 


6713.74 


-2.250e+02 


3.797e+01 


1.056e-01 


1.325e+01 


-1.090e-05 


-3.997e-05 


4647 


6161 


-0.32 


0.39 


133 


0.019 


193 


6725.36 


-2.724e+()2 


3.980e+01 


1.287e-01 


7.781e+00 


-1.362e-05 


9.816e-05 


4698 


5989 


-0.26 


0.39 


134 


0.020 


382 


6726.67 


-1.289e+01 


4.785e+()l 


4.599e-02 


2.696e+00 


-6.151e-06 


5.169e-04 


4671 


6403 


-0.78 


0.39 


135 


0.017 


311 


6733.15 


-2.182e+02 


4.476e+01 


1.088e-01 


7.913e+00 


-1.148e-05 


2.804e-05 


4624 


6227 


-0.45 


0.39 


136 


0.019 


172 


6739.52 


2.370e+01 


6.136e+01 


3.417e-02 


9.389e+00 


-6.356e-06 


-4.288e-03 


4556 


5827 


-0.72 


0.39 


137 


0.021 


306 


6786.86 


-1.150e+02 


5.132e+01 


7.364e-02 


8.126e+00 


-8.511e-06 


-1.370e-03 


4556 


6276 


-0.54 


0.39 


138 


0.024 


110 


6793.26 


4.719e+01 


2.173e+02 


5.951e-04 


5.214e+01 


-1.032e-06 


-3.675e-02 


4647 


5914 


-0.20 


0.39 


139 


0.027 


359 


6806.85 


-4.387e+()l 


4.316e+01 


6.488e-02 


4.865e+00 


-8.845e-06 


9.710e-04 


4556 


6361 


-0.78 


0.35 


140 


0.021 


376 


6810.27 


-4.150e+01 


3.646e+()l 


5.659e-02 


1.603e+00 


-7.038e-06 


2.640e-03 


4595 


6374 


-0.84 


0.35 


141 


0.016 


377 


6820.37 


-2.306e+01 


5.392e+01 


4.775e-02 


7.352e+00 


-6.336e-06 


-3.263e-04 


4649 


6374 


-0.78 


0.35 


142 


0.027 


367 


6828.60 


-1.282e+02 


1.418e+01 


8.794e-02 


1.754e+00 


-9.707e-06 


6.572e-03 


4556 


6403 


-0.84 


0.33 


143 


0.024 


174 


6837.01 


-1.395e+02 


3.806e+()l 


7.098e-02 


1.724e+01 


-7.564e-06 


-1.265e-03 


4556 


6161 


-0.22 


0.39 


144 


0.028 


345 


6839.84 


2.072e+02 


1.071e+02 


-2.714e-02 


1.279e+01 


-5.973e-07 


-1.090e-02 


4556 


6287 


-0.78 


0.39 


145 


0.019 


384 


6842.69 


-7.707e+01 


6.390e+01 


6.450e-02 


2.873e+00 


-7.684e-06 


-2.778e-03 


4671 


6289 


-0.70 


0.35 


146 


0.025 


328 


6843.66 


-2.222e+02 


3.112e-01 


1.244e-01 


-1.969e+00 


-1.305e-05 


9.038e-03 


4595 


6374 


-0.84 


0.33 


147 


0.033 


205 


6855.72 


1.443e+02 


2.005e+02 


-3.124e-02 


5.407e+01 


1.691e-06 


-3.082e-02 


4595 


6161 


-0.25 


0.39 


148 


0.017 


272 


6857.25 


-1.797e+02 


1.675e+01 


9.315e-02 


8.259e+00 


-1.005e-05 


4.228e-03 


4556 


6227 


-0.52 


0.39 


149 


0.028 


379 


6858.15 


-1.407e+02 


6.900e+00 


9.054e-02 


-6.242e-()l 


-9.878e-06 


7.621e-03 


4556 


6374 


-0.84 


0.35 
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Table 2. continued. 



Sousa, S. G. et al.: TMCalc - a fast code to derive Teff and [Fe/H] for FGK stars 



id 


stddev 


ns 


A 


co 


c\ 


C2 


C3 


c 4 


C5 


eff 


rrmax 

eff 


fehmin 


fehmax t 


150 


0.025 


244 


6861.94 


2.667e+02 


9.333e+01 


-5.126e-02 


2.030e+01 


1.418e-06 


-8.969e-03 


4578 


6026 


-0.72 


0.39 


151 


0.026 


340 


6862.50 


-2.833e+02 


2.371e+01 


1.330e-01 


7.602e+00 


-1.359e-05 


3.926e-03 


4556 


6361 


-0.54 


0.39 



t The columns correspond to the following: id is the ID number of the line; ns is the number of stars used for line calibration (out of a total of 451); stddev is 
the final standard deviation for each calibration; A is the wavelength of the calibrated line; Co, C1.c2.c3, C4, and c; are the fitting coefficients for each line; 7/™" and 
T™fi are the limits in effective temperature of the stars used to calibrate each line; fehmin and fehmax are the limits in [Fe/H] of the stars used to calibrate each line. 



10 



